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The interaction between a wide ultrarelativistic fully-ionized plasma beam and a magnetic barrier 
is studied numerically. It is assumed that the plasma beam is initially homogeneous and impacts 
with the Lorentz factor Fq 3> 1 on the barrier. The magnetic field of the barrier Bq is uniform and 
transverse to the beam velocity. When the energy densities of the beam and the magnetic field are 
comparable, a — ^unompt? (Fo — 1)/ Bq ~ 1, the process of the beam - barrier interaction is strongly 
nonstationary, and the density of reversed protons is modulated in space by a factor of 10 or so. The 
modulation of reversed protons decreases with decrease of a. The beam is found to penetrate deep 
into the barrier provided that a > ttcr, where Ocr is about 0.4. The speed of such a penetration is 
subrelativistic and depends on a. Strong electric fields are generated near the front of the barrier, 
and electrons are accelerated in these fields up to the mean energy of protons, i.e. up to ~ mpC?Vo. 
The synchrotron radiation of high-energy electrons from the front vicinity is calculated. Stationary 
solutions for the beam - barrier collision are considered. It is shown that such a solution may be 
only at a < 0.2 — 0.5 depending on the boundary conditions for the electric field in the region of the 
beam - barrier interaction. Some astrophysical applications of these results are briefly discussed. 
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I. INTRODUCTION 

There is now compelling evidence that plasma is 
ejected from many astronomical objects and flows away 
at relativistic speeds. The Lorentz factor, Fq, of such 
plasma wind is a few x(l — 10) for the jets associated 
with active galactic nuclei [P and ^ 10^ — lO'^ or even 
more for the 7-ray bursters B . A strong magnetic field 
may be in the outflowing gas [^|-^. Relativistic magne- 
tized winds can interact with an external medium (e.g., 
an ordinary interstellar medium). It was pointed out 
in Ref. Q that such an interaction may be responsible 
for radiation of both X-rays and 7-rays from the 7-ray 
bursters. For consideration of the interaction between a 
relativistic magnetized wind and an external medium, it 
is convenient to switch to the wind frame. In this frame, 
the problem of the magnetized wind - external medium 
interaction is identical to the problem of collision between 
a wide relativistic beam of cold plasma and a region with 
a strong magnetic field which is called a magnetic barrier. 

The problem of the interaction between plasma beams 
and magnetic barriers was attacked in many experi- 
ments and theoretical papers (see Ref. |^ and references 
therein). However, all known experimental studies of the 
beam - magnetic barrier interaction are irrelevant to our 
problem. This is because the plasma beams produced by 
the laboratory equipment are either non-relativistic or 
very narrow and the cross-beam sizes are much smaller 
than the gyroradius of the beam protons in the field of 
the barrier. In the theoretical studies, it is usually used 
the assumption of Rosenbluth that the charge separa- 
tion is small. However, this assumption is valid provided 
that the inequality (wq/c)^ ^ rrie/mp ~ 5.4 x lO^"* holds 



true, where vq is the beam velocity, me is the electron 
mass and rup is the proton mass. For relativistic beams, 
vo ~ c-\/l — (1/Fo)^^ ~ c, the assumption of Rosen- 
bluth proves inadequate, and the charge separation is 
very important for the dynamics of particles near the 
barrier front pO|jTl| ]. Below, the interaction of a wide 
ultrarelativistic plasma beam with a magnetic barrier is 
considered numerically (for preliminary results and some 
astrophysical applications, see Ref. pd|). 



II. FORMULATION OF THE PROBLEM AND 
BASIC EQUATIONS 



The situation to be discussed is the following. At the 
initial moment, t = 0, the ultrarelativistic (Lorentz factor 
Fg ^ 1) neutral beam of protons and electrons (number 
densities Hp = = hq) runs along the x axis into the 
magnetic barrier which is the half-space x > with an 
external magnetic field Bq = BQezQ[x], where uq, Bq are 
constants and Q [x] is the step function equal to unity for 
a; > and to zero for x < 0. The beam is infinite in the y 
- z dimensions and semi-infinite in the x dimension. Our 
goal is to construct a l^D time-dependent solution for 
the problem, i.e. to find induced electromagnetic fields 
(E = E^e^ + EyBy-, B = Bbz) and motion of the beam 
particles in the x ~ y plane. The field structure and the 
beam particle motion are to be treated self-consistently. 
All the quantities are assumed to be dependent on t and 
X only. 

As noted earlier, the strength of the magnetic fields in 
the astrophysical winds may be very high, especially in 
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the winds outflowing from the 7-ray bursters |4||^ . High- 
energy electrons generate synchrotron radiation in these 
fields, and the radiation damping force which acts on the 
radiating electrons has to be taken into account. Since 
nip ^ TOe, both the synchrotron radiation of protons and 
the radiation damping force which acts on the protons are 



very small (e.g., Ref. |12|) and may be neglected. 

The following set of equations can be used to describe 
the process of the beam - barrier collision II4I : 
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where p and j are the densities of charges and cur- 
rents respectively, e is the charge of particles, m is the 
mass of particles, m = {mg, mp}, s is the interval, 
— {r, (v/c)r} is the four- velocity, F is the Lorentz 
factor of particles, and F"^^ is the electromagnetic field 
tensor. No simplifying assumptions besides geometrical 
ones are exploited when rewriting Maxwell equations (|l|) 
- (^). The second term in the right-hand side of the rel- 
ativistic equation of particle motio n (|4| ) is the radiation 
damping four-force. Following Ref. only the term of 
the highest power in F is left in the damping force. This 
is valid for ultrarelativistic motion of particles, F 3> 1 
(e.g., Ref. dl). 

The four-velocity from Eq. (^ is related to the 
charge and current densities which appear in the Maxwell 
equations ([T]) - (^: 



f = {cp,3x,3y.^} 



(5) 



where ^ means the sum over the beam particles per uni- 
tary volume. 

The energy losses, /, of electrons because of their ra- 
diation in the electric and magnetic fields are computed 
using the following equation p2| : 



/ = 




imic 



(6) 



To evaluate the spectrum of the radiation of electrons, 
we have used the following expression for the spectral 
intensity of synchrotron radiation M]: 



lu = 2 / K5/3{r])dr], 



(7) 



where i^s/s is the modified Bessel functions of 5/3 order. 
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(8) 



In Eqs. (^ and (||), it is taken into account that in 
our case the velocities of particles are perpendicular to 
the magnetic field, v _L B. 

The initial beam number density is taken from the 
equation: 



nompC (To - 1) = a — 
on 



(9) 



where a is the dimensionless parameter. In this work 
we study mainly the case a ^ 1 when the plasma fiow 
pressure is comparable with the magnetic field pressure 
of the barrier. In different runs, a ranges from 0.2 to 4. 

To integrate the set of equations (||) - (||), we used 
a macro-particle approximation in which all particles of 
the beam are subdivided into a large number, ^ 10'', 
of packets. The particles within a packet are bundled 
together forming a large "macro-particle" which is infini- 
tively small along the axis x and infinitively large in all 
directions of the y, z plane. Each packet contains either 
protons or electrons (for details of the numerical method, 
see Ref. |l|). 

The examined space-time domain is the following: 



< X < x„ 



and <t <t„ 



few x(l - 10)(c/wBp), Xn 



(10) 



a few 



where Xmin = 

x(l-10)(c/wsp), imax = a fcw x(l-10)Tp, Tp = 2Tr/u;Bp 
is the proton gyroperiod, and LUBp = eBo/{mpCTo) and 
c/ojBp = "ipC^Fo/(ei3o) are the proton gyrofrequency 
and gyroradius, respectively. The time step of calcula- 
tions is a few xlO^^Tp, depending on a. 

The boundary condition for E^ is E^ |a:=2;max = 0, as 
a^max is chosen so that no beam particles penetrate deeper 
than Xniax in the observed interval of time. Electromag- 
netic waves generated due to charge flow are allowed to 
escape freely from the system that is used as boundary 
conditions for B and Ey in Eqs. (||) and (||). 



III. RESULTS OF NUMERICAL SIMULATIONS 
AND SCALING 



The relativistic cold plasma - magnetic barrier colli- 
sion is characterized by the following parameters: Bq, 
Tq and a. The input parameters of our simulations are 
given in Table ||. The values of Bq and Fq are chosen to 
be relevant to cosmological 7-ray bursters P-p|j7|,pT[| . 
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A. Particle dynamics and penetration of the beam 
particles into the barrier 

The density of protons which move towards the mag- 
netic barrier, Vx > 0, is almost unperturbed untill the 
distance to the barrier is smaller than ^ c/usp- Re- 
versed protons which move away from the magnetic bar- 
rier, Vx < 0, are bunched in the process of the beam 
- barrier interaction (Fig. The modulation of the 
density of reversed protons increases with increase of a. 
Namely, the ratio of the maximum to minimum densities 
of reversed protons is '--^ 2 at a ~ 0.2 and ~ 10 at a ~ 1. 
A typical length of such a modulation is roughly the pro- 
ton gyroradius c/ujBp- Similar phenomena was observed 
also in numerical simulations of coUisionlcss shock waves 
near the shock front (e.g., Ref. @). 

At a ~ 1, the mean Lorentz factor of reversed protons 
that far enough from the barrier, x < —(c/ujBp), where 
the process of strong interaction between particles and 
fields is more or less over, is (r°"*) ~ (0.5 ± 0.1)ro, i.e. 
about half of the initial kinetic energy of protons is lost in 
the process of their collision with the barrier (see Table 
U) . Figure || shows the energy distribution for reversed 
protons. 

The density of electrons which move away from the 
barrier is well correlated with the density of reversed 
protons; the correlation coefficient is r ~ 0.8, where 
r = o-pe/y/cTppaee, and 



-){nj-nj)dx, {i,j} = {e,p}. (11) 



These electrons screen the electrostatic field. Ex, of pro- 
ton bunches mostly, but not completely. 

At a > Qfcr, ttcr — 0.4, it is observed that the length of 
the beam particle penetration into the barrier increases in 
time. Figure || shows the x coordinate, Xpcn, of the most 
deeply penetrated proton as a function of time t. The 
speed of such a proton along the x axis varies strongly one 
from another. However, the mean velocity, Vpcn = i^x), 
of the penetrated protons remains more or less constant 
within the studied time intervals, where (vx) is the pro- 
ton velocity which is averaged over proton gyroperiod. 
The value of Vpon depends on a and is equal to zero (no 
penetration) at a < aa — 0.4 (see Fig. ^ and Table |l|). 

In the region x < Xpcn, strong longitudinal (Ex) and 
transverse {Ey, Bz) electromagnetic waves propagate 
while almost no remnants of the external magnetic field 
BoQN are found (Fig. ^). Roughly, it could be stated 
that at a > ctcr the magnetic barrier is pushed according 
to the law B{t) = BqQ[x — Xpc-ci{t)\ and a;pon appears to 
be a location of the barrier front at the moment t as well 
as the particle penetration depth. At x > Xpcm the mag- 
netic field remains unchanged except for its time-space 
variations due to low-frequency electromagnetic waves 
which are generated by the time-variable currents j in 



the front vicinity. At a ~ 1, the typical amplitude of 
these waves is Bq ~ (0.2 — 0.3)i3o- 

In the case a < ctcr, when the ion penetration length 
does not increase with time, all protons move along the 
same track with small deviations from it. In this case, 
the velocity of all reversed protons is perpendicular to 
the barrier front. However, at a > acr the trajectories 
of protons differ qualitatively one from another (for the 
angular distribution of reversed protons in this case, see 
Fig. H). Therefore, the dimensionless density a may be 
called a "stochastization parameter" of the beam pro- 
tons. In all runs, the trajectories of electrons are very 
chaotic and differ one from another qualitatively. 



B. Acceleration of electrons and high-frequency 
radiation 

The most important feature of the dynamics of elec- 
trons is that they are accelerated in the barrier front 
vicinity by induced electric fields and, thus, accumulate 
substantial portion of the kinetic energy of protons (see 
Fig. I and Table p. Figure M shows the energy spec- 
trum of outflowing electrons, Vx < 0, in run N. At a ^ 1, 
the mean Lorentz factor of outflowing electrons and their 
maximum Lorentz factor far enough from the barrier, 
X < —(c/uBp), are 
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-0.05 (^) To, r:-,,c.lp)ro (12) 



within a factor of 2. The fraction of the kinetic energy of 
relativistic protons which is transformed into the energy 
of outflowing electrons is up to ~ 10%. The rest of the en- 
ergy that is lost by protons is transformed into both low- 
frequency electromagnetic waves and synchrotron high- 
frequency radiation. 

The synchrotron high-frequency radiation is generated 
by single electrons in the process of their ultrarelativis- 
tic motion in the electromagnetic fields. The mean en- 
ergy of so-called radiating electrons near by the barrier 
front, X > ~{c/lubp), is several times higher than the 
mean energy of outflowing electrons far from the barrier, 
X < -{c/ujBp), i.e. {Tl^'^) ~ (3 - 5)(r™*) (see Table 
II). The radiating electrons are responsible for genera- 
tion of the main part of the synchrotron high-frequency 
emission. 

If the average fraction, of the kinetic energy of the 
plasma beam that is radiated in the vicinity of the mag- 
netic barrier is small, < 0.01 or TqBq < 3 x 10® G, the 
results of our simulations for a ~ 1 may be fitted by the 
following analytic expression: 



^^~3x 10-4(ro/102)2(Bo/103G). 



(13) 



At LgSo ^ 10® G, the value of tends asymptotically 
to ^ 0.15 as FqSo increases. The characteristic energy 
of synchrotron photons is 
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(e^) ~ 3(^^/10-2) MeV 



(14) 



within a factor of 2 or so. Except for numerical factors, 
Eqs. (|l3| ) and ( p^ can be derived analytically from Eqs. 
(i-i)- 

At a <C 1, electromagnetic fields which are induced 
in the process of the beam - barrier collision are much 
smaller than Bq, and the beam particles may be treated 
as test particles in the field of the barrier. In this case, 
there is almost no energy transfer from protons to elec- 
trons. Therefore, it is natural that both acceleration of 
electrons and their high-frequency radiation are strongly 
suppressed at a <C 1 as was observed in our simulations 
(see Table Preliminary, we have the following scaling 
laws: (r°"*) cx a, r°^;^^^^ cx a, (r'^^^) cx a, cx and 
(e-y) oc a^. 

We did not observe any decrease in both non- 
stationarity of the process of the beam - barrier interac- 
tion and in acceleration of electrons, even in simulations 
extended up to tmax — 16Tp (see Figs. |l| and ||). 



IV. STATIONARY COLLISION 

The problem of the beam - barrier collision is sim- 
plified significantly if it is treated as a stationary one. 
This is because in a stationary consideration a time- 
dependence of all values is disregarded, and the beam 
particles of any kind (protons or electrons) move along 
exactly the same track. In this Section we consider the 
beam - barrier interaction in a stationary manner. 

A fully relativistic self-consistent model of stationary 
interaction between a wide plasma beam and a magnetic 
barrier was developed by Peter, Ron and Rostocker in 
their pioneering paper |10[| . The set of equations which 
describes such an interaction is m(i 
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rrie/mp 



1 



dE^ 
dx 



— —Err, 



dA 



'-^B^. 



dx ' dx 

where F^ = Fq - tijcf), z.j = F^ - ^I'^-A^ - 1, 



Mi 



— 1 for electrons (j = e) , 

TOe / mp for protons (j = p) , 



(15) 



(16) 



(17) 



(18) 



(f> and Ay are the electrostatic and magnetic potentials 
in units of e/mc^. In these equations we have mea- 
sured the distance x in terms of \/2('^pe/3o/c)) where 
LOpe = (47rnoe^/me)^^^ is the electron plasma frequency 
and /3o = fo/c 



Eqs. dH) and (|T|) for the magnetic and electric fields 
are first order ordinary differential equations. Strictly 
speaking, such an equation allows only one boundary con- 
dition. For the magnetic field B^ , the boundary condition 
is Bz = Bq at the proton turning point Xp (see Fig. ||), 
where Bq is the field of the barrier. The strength of the 
magnetic field Bf at the front of the barrier, x = 0, has 
to be found from integration of Eq. (15). An important 
issue is where the electric field E^ should vanish: at the 
right or left boundary of the collision regin < x < Xp. 
In the paper [|lO|, the electric field vanishes at x = Xp, 
Ex\x=xp = 0. This seems quite natural and corresponds 
to our nonstationary solution. However, in Ref. |10| it 
is wrongly suggested that if the beam density is high 
enough the electric field E^ is zero at the left boundary 
as well, Ex\x=o = 0. Below, stationary solutions for the 
beam - barrier collision are considered in detail, and, in 
particular, it is shown that there is a net charge in the 
region < a; < Xp irrespective of the beam density un- 
till the beam - barrier collision is stationary, i.e. the 
transverse electric field Ex cannot be zero at the both 
boundaries at once. 

The right-hand side of Eqs. (|l5|) and ( [l6| ) for the 
derivatives of the electric an magnetic fields contains 
terms which tend to infinity at the electron and proton 
turning points, i.e. at both x = x,, and x = Xp. This 
is the main difficulty which does not permit to apply a 
standard numerical technique to integration of Eqs. ( p^ ) 
and (^6|) directly. Since the electric and magnetic fields 
should remain finite everywhere, the mentioned singular- 
ities are integrable. To integrate a set of equations with 
such a weak, integrable singularity, one should switch to 
another independent variable so that all derivatives with 
respect to this new variable are finite. The interval of 
the beam particles may be used as such a new indepen- 
dent variable instead of x. In this case, the right-hand 
sides of Eqs. (|l5|) and (^6|) are finite everywhere. This is 
because the particle density per unitary interval is con- 
stant (e.g., Ref. ||l^), rather than infinite at the turning 
points if measured per unitary x. We integrated Eqs. 
( |l5| ) - (0) from x — Qtox — Xe using the electron inter- 
val Sg as an independent variable. Then, we switched to 
another independent variable, Sp, which is the proton in- 
terval measured from the electron turning point with the 
initial condition Sp\x=x^ = Se\x=x^- We thus broke the 
region of integration into two parts and used a different 
independent variable for each sub-interval to integrate 
Eqs. (|l^) - (|l7|). Then, we switched back to x as an in- 
dependent variable to present results of this integration. 

Following Ref. [nO|, we have used the dimensionless 



l/B? as a measure of 



parameter ctf — STrno^^pC^yFj — ^/jj^ 
the beam density instead of a. The parameter a/ is more 
convenient than a for integration of Eqs. (15) - (17) be- 
cause such integration is performed from x = Qto x = Xp. 
In our calculations the value of varies in a wide range 
from 10~2 to 10^. 
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To integrate Eq. ([l6|), the following two kinds of a 
boundary condition for were specified: E^I^^q = or 
^x\x=xp — ^- Integration of Eq. for the boundary 
condition E^l^^^ = was carried out with a "shoot- 
ing method" technique. The electric field was guessed 
at the left boundary, x — 0, and Eq. ( [l6| ) was inte- 
grated with all the quantities specified on that bound- 
ary. Then, the electric field on the left boundary was 
adjusted to give a closer agreement with the boundary 
condition Ej;\^^^ — 0. Such shots were repeated untill 
the boundary condition for the electric field was satisfied 
with a desirable accuracy. As a rule, only a few shots 
were required to reach the accuracy of ^ 1%. 

The computational engine we used to solve the regular- 
ized set of equations was "LSODE" package [|l5| together 
with the "OCTAVE" interpreter The developed oc- 
tave script is mostly compatible with the popular MatLab 
language. 

Figure ^ shows typical numerical solutions for the lon- 
gitudinal electric field E^, the magnetic field and the 
X component of the four-velocity of electrons Ue^x iu run 
SRIO. In this run, the boundary condition for E^ is 
Ex\x=xp = 0. The energies of both electrons and pro- 
tons inside of the barrier are shown in Fig. ^ The 
results of our calculations are summarized in Tables IV 
and V. From these Tables and Fig. ^ we can see that 
the electric field Ex is not equal to zero at x = and 
X = Xp at once, as mentioned above. A common fea- 
ture of all stationary solutions is that the electric and 
magnetic fields sharply increase near the electron turn- 
ing point Xe as it was pointed out in Ref. [0. There is 
the following tendency: the greater the value of a/ is, 
the sharper the rise of the fields is. Due to the change of 
independent variables discussed above, the sharp rise of 
the fields at the electron turning point is well-resolved: 
it covers a large number of integration points and deriva- 
tives of the fields with respect to the new independent 
variable remain finite (see Fig. ^l|). The later facts allow 
us to claim a high accuracy of our calculations. 

In the paper of Peter et al. |lo), Eqs. ((ll) - (|l|) 
were integrated directly using x as an independent vari- 
able, and, therefore, they ran into serious difficulties on 
their way which led to certain mistakes. For exam- 
ple, they were not able to calculate the value Exjx) at 
X = Xe directly, and the integration of Eq. (Hq) was 
treated as a two-point boundary value problem, adjust- 
ing Ex{xe) to be equal to zero at the both boundaries, 
-£'k(0) = Ex{xp) = 0. Our calculations show that such 
an adjustment is impossible (see Tables IV and V). The 
authors of Ref. ||f^ believed that at a/ ^ 1 the orbit of 
electrons is strongly elongated in the y direction in a nar- 
row vicinity of the electron turning point, and the beam 
electrons move more or less along the barrier front with 
X ~ Xe for a long time. As a result, a net charge in the 
collision region might be zero. We did not observe such a 



strong elongation of the electron orbit. Moreover, for the 
boundary condition Ex\^^q = we have a/ ^ 1, i.e. for 
this boundary condition a/ cannot be much more than 
unit. 

We have confirmed the suggestion of Ref. that elec- 
trons of the beam may be strongly accelerated inside the 
barrier. However, such an acceleration was observed only 
in runs with the boundary condition for Ex in the form 
Ex\^^^ = 0. For the boundary condition -E^I^^q = 0, 
the beam electrons do not penetrate substantially deeper 
into the barrier than their gyroradius in the field Bf, and 
their acceleration is very small. 

For the boundary condition Ex\x=xp = 0, a stationary 
solution of Eqs. (|l^) - ( |l7| ) can exist for any nonzero 
value of oj. For all these solutions, we have a 0.5 
(see Fig. |2|). This upper limit on a is more or less ev- 
ident since in a stationary case the ram pressure of the 
beam cannot be more than the magnetic field pressure 
of the barrier. As to the case of the boundary condi- 
tion Ex\x=Q = 0, a stationary solution can exist only if 
ck/ 1 and a < 0.2. This reduction of the upper limit 
on a occurs because in such a solution the electric field 
Ex far enough from the barrier front is directed inside 
the barrier, E'^, > at x ^ a few xxg. This field causes 
protons to penetrate even deeper into the barrier and to 
generate even stronger electric field. At a/ > 1, at some 
distance from the barrier front the electric field inside the 
collision region is stronger than the magnetic field, and 
the barrier cannot suppress penetration of protons inside 
itself. 



V. DISCUSSION 

One of the main results of our time-dependent simu- 
lations which are presented in Sec. Ill is that the out- 
flowing electrons are strongly accelerated near the barrier 
front and accumulate substantial portion of the kinetic 
energy of protons. At a ^ 1, when the ram pressure 
of the beam is equal about the magnetic field pressure 
of the barrier, the fraction of the kinetic energy of pro- 
tons that is transferred to outgoing electrons is up to ~ 
10%. The acceleration of outflowing electrons is com- 
pletely due to nonstationarity of electromagnetic fields 
which are induced inside and near the barrier during the 
beam - barrier interaction. Indeed, in a stationary case 
the electric field is constant in time and depends on the x 
coordinate only, E = {Ex{x),{), 0}. If, as in both Ref. 0] 
and Sec. IV, the energy losses of electrons due to their 
radiation inside the barrier are ignored, the electron en- 
ergy is a function of x only, and the energy of reversed 
electrons outside the barrier coincides with their initial 
energy before the beam - barrier collision. 

The bunching of reversed protons is a key element of 
nonstationarity of the beam - barrier collision. Roughly, 
the process of bunch formation may be illustrated in the 
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following way. The first protons, which are entered into 
the barrier from t = to t <C Tp, run along almost semi- 
circle trajectory in almost unperturbed field of the mag- 
netic barrier. For the same Lorentz-factor, the gyrora- 
dius of electrons is much smaller than the gyroradius of 
protons. Therefore, electrons can not penetrate as deep 
as protons into the barrier. The separation of electric 
charges induce a strong electric field in the x direc- 
tion according to Eq. (|l|). The following protons feel 
this electric field. The field Ex decelerates protons and 
accelerates electrons. Protons that are injected into the 
barrier later run along "shorter" trajectories and spend 
less time inside the barrier. As a result, in about half of 
a proton gyroperiod, t ~ ^Tp, most of the protons quit 
the barrier almost simultaneously, and formation of the 
second bunch begins, and so on. 

In case of large a, the velocities, Vx, of reversed pro- 
tons in the x direction are quite different (see Fig. ||) 
so that the bunches will decay at some distance from the 
barrier. Once there is no bunching at all in low a case, no 
bunches can be observed far from the barrier no matter 
what a is considered. 

The propagation of plasma beams across a magnetic 
field is one of the oldest problems in plasma physics. In 
spite of the long history of investigation (e.g., Ref. [^), a 
clear model of this phenomenon has yet to emerge. How- 
ever, some general conclusions about the beam - mag- 
netic field interaction were done many years ago. For 
example, if the plasma flow pressure is more than the 
magnetic field pressure, it is possible for the plasma to 
cross the field by purely magnetohydrodynamic princi- 
ples. This propagation mode is equivalent to the motion 
of a solid conductor across the field. During propagation, 
the plasma beam picks up and carries along the ambient 
plasma and magnetic field. Such models of the beam 
propagation into the magnetic field are based on the dia- 
magnetic properties of the plasma. Most probably, in our 
time-dependent simulations we observed such a propaga- 
tion mode. Physically, the pushing of the magnetic field 
is provided by strong transverse electric currents (jy) 
in the domain Xpcn(i) — (c/ujBp) < x < Xpon(t), where 
Xpcn{t) is the X coordinate of the most deeply penetrated 
particle. If a is large enough, a > a^, this current layer 
screens strongly the external magnetic field of the barrier 
at a; < Xpon(t) — (c/uBp), and moves deeper and deeper 
into the barrier. From our time-dependent simulations, 
we have aa- ~ 0.4. This is more or less consistent with 
the conclusion that there are no stationary solutions for 
the beam - barrier collision if a is more than ^ (0.2 — 0.5) 
depending on the boundary condition on the field Ex- 

Many important features of our nonstationary solu- 
tions are absent in a stationary one. They are energy 
transfer from incoming protons to outgoing electrons, 
bunching of outgoing protons, excitation of low- frequency 
electromagnetic waves, etc. 

The main question which remains open untill now is 



what the final state of the beam - barrier system in a very 
long-time run, t ^ Tp, is. As noted, we did not observe 
any decay of non-stationarity of the beam - barrier inter- 
action, even in simulations extended up to tmax — IGTp. 
However, it is possible that at t ^ lOTp our nonstation- 
ary solution tends to a stationary one. It is worth noting 
that a stationary solution of the beam - barrier collision 
might be unstable if a is high enough. Indeed, the en- 
ergy of electrons inside the barrier may be a few hundred 
times higher than their initial energy. In this case, small 
(~ 10~^) perturbations might result in that some part 
of reversed, Vx < 0, electrons is stopped and captured 
inside the barrier before reaching the barrier front. Be- 
sides, electrons may be captured by the barrier because 
of their energy losses via synchrotron emission in electric 
and magnetic fields if these fields are strong enough. In 
turn, such a capture of electrons may result in complete 
destruction of the stationary solution. 

The beam - barrier collision is similar, in many re- 
spects, to collisionless shocks. Such shocks in astrophys- 
ical settings can and do accelerate charged particles to 
high energies (see, for a review, Refs. [ p^pc| ]). The effi- 
ciency of particle acceleration by shocks may be as high 
as - 20 - 40 % (e.g., Ref. ||l|l). This value is more or 
less the same as the fraction of the kinetic energy of the 
beam which may be transferred to high-energy electrons 
in the process of the beam - barrier collision. 

Besides of analytical calculations (e.g., Refs. | |2l| , p^ ), 
acceleration of particles by shocks was studied numer- 
ically in both test particle Monte Carlo simulations 



(see Ref. |23 and references therein) and self-consistent 
plasma simulations (e.g., Refs. p4|-[27|]). Most current 
self-consistent simulations of plasma shocks used in astro- 
physics are of the hybrid type because the rrie/mp ratio 
is small, me/mp ~ 0.54 x 10"'^. In the hybrid approach, 
the ions are treated kinetically using standard particle- 
in-cell techniques, while the electrons are treated as a 
massless, charge neutralizing fluid. In our simulations 
we did not use the hybrid approach because the case 
of ultrarelativistic beam is, in some respects, easy than 
non-relativistic case. This is because in our case, Fq 3> 1, 
electrons are accelerated fast, and the mean energy be- 
comes only an order of magnitude smaller than the mean 
energy of protons. In this case, the ratio of the gyroradii 
of electrons and protons is about two orders of magni- 
tude larger than rrie/mp that is the ratio of the gyroradii 
of electrons and protons in non-relativistic plasma. This 
allowed us to treat both electrons and protons kinetically. 

The configuration considered in our simulations is not 
quite self-consistent. Indeed, the barrier magnetic field 
has to be generated by some kind of current flowing in 
the plane x — in the y direction. However, the evo- 
lution of this current because of its interaction with the 
induced electromagnetic fields was not studied. In labo- 
ratory experiments, an external current flowing along a 
sheet of thin wires can be the barrier-front current. As to 
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a situation which is relevant for astrophysics, the barrier 
fields has to be generated by the magnetized plasma of 
the barrier. Numerical considerations of such a configu- 
ration are under way. 

VI. ASTROPHYSICAL APPLICATIONS 

There are at least two kinds of astrophysical objects 
to which the results of our simulations may be applied. 
They are the sources of 7-ray bursts and strongly coUi- 
mated jets of active galactic nuclei. 

A. Cosmological 7-ray bursters 

Gamma-ray bursts are brief, ~ 10^^ — 10^ s, bursts 
of high-energy radiation that appear at random in the 
sky, emitting the bulk of their energy in the range from 
~ 0.1 MeV to a few MeV (for a review of observational 
data on 7-ray bursts see Ref. ||2^ ) . During their appear- 
ance, they often outshine all other sources in the 7-ray 
sky combined. More than two decades have passed since 
the discovery of 7-ray bursts, but their origin is still a 
mystery. The observed isotropy of 7-ray bursts in the sky 
and deficiency of faint bursts strongly suggest that 
they are cosmological in origin ||3l[] . A cosmological origin 
of 7-ray bursters implies that the typical energy output 
in both hard X-rays and 7-rays is about a few xlO^^ 
ergs, that is, about IO'^Mqc^. Such high energetics of 
7-ray bursters and a short time scale of 7-ray flux vari- 
ability call for very compact objects as a source of 7-ray 
bursts. These objects may be either fast-rotating neu- 
tron stars which arose from accretion-induced collapse of 
white dwarfs in binaries or differentially rotating disk- 
like objects which are formed by the merger of a binary 
consisting of two neutron stars ||3^,Q. The strength of 
the magnetic field Bg at the surface of these objects may 
be up to - IQis G |,|. 

A common point of all acceptable models of 7-ray 
bursters is that a relativistic, Tq > 10^ — 10"^, wind is a 
source of their radiation (e.g., Refs. [t34y35[ ). Otherwise, 
high-energy photons are absorbed by 7 -I- 7 ^ e+ -I- e~, 
and both observed durations of 7-ray bursts and their 
hard energy spectra, often with a very high energy tail 
extending up to hundreds of MeV, cannot be explained. 
Most probably, the radiation of 7-ray bursts is produced 
in the process of interaction between such a wind and an 
external medium 0Jll[l. For both the post-merger objects 
and fast-rotating neutron stars, the total energy of the 
outflowing gas cannot be more than 10'^'^ ergs. Therefore, 
to explain the energetics of cosmological 7-ray bursters 
it is needed that not less than a few per cent of the wind 
energy is transformed into high-frequency radiation. 

From our simulations of the beam - magnetic barrier 
collision it follows that the energy which is lost by a rela- 



tivistic magnetized wind in the process of its interaction 
with an external medium is distributed in the following 
way: About half of this energy is in ultra-relativistic pro- 
tons which are reflected from the wind front. The mean 
energy of reflected protons in the frame of the burster is 
~ iripC^TQ. The other half of the wind energy losses is 
distributed more or less evenly between low-frequency os- 
cillations of electromagnetic fields, and both high-energy 
electrons and their high-frequency radiation. The low- 
frequency electromagnetic oscillations which are gener- 
ated near the wind front can accelerate electrons to very 
high energies and decay at some distances from the front 
[^,0. As a result of these, the fraction of the kinetic en- 
ergy of relativistic magnetized winds which is transferred 
to high-energy electrons and then to high-frequency emis- 
sion might be, in principle, as high as ~ 40 — 50%. This 
is more than enough to explain the energetics of 7-ray 
bursters. 

In the burster frame, the characteristic energy, (e-y), of 
photons which are generated by accelerated electrons in 
the vicinity of the wind front is ^ Fq times higher than 
the value of (e^) (see Eq. (p^) in the frame of the out- 
fiowing gas because of the Dopplcr effect, (e^) ~ ro{e.y). 
Taking into account that Fq has to be more than ^ 10^ 
[^,||, from Eq. ( p^ we can see that the radiation which 
is generated near the wind front with a high enough ef- 
ficiency, ^.y ^ 10~^, has the characteristic energy {e.y) ^ 
10^ — 10'^ MeV, while the main part of emission of de- 
tected bursts is in the range from a few x 10 keV to a few 
MeV. In other words, the radiation from the front vicin- 
ity is either too hard or too weak to explain the bulk of 
the burst emission. The 7-ray emission which is gener- 
ated near the wind front may be responsible for the high 
energy 7-rays, e^y > a few ten MeV, which are observed 
in the spectra of some bursts 

The bulk of emission from 7-ray bursters may be ex- 
plain by the synchrotron radiation of high-energy elec- 
trons which are ejected from the wind front ahead of 
the front. Indeed, electromagnetic fields ahead of the 
wind front are the fields of low-frequency electromagnetic 
waves which are generated near the wind front due to 
non-stationarity of the wind ~ external medium interac- 
tion. In this case, the mean value of the magnetic field is 
about 0.2_Bo. The mean energy of high-energy electrons 
reflected from the wind front is ~ (3 — 5) times smaller 
than the mean energy of radiating electrons which are 
responsible for the generation of 7-rays near the front. 
Therefore, the characteristic energy of synchrotron pho- 
tons which are generated far ahead of the wind front is 
about hundred times smaller than the characteristic en- 
ergy, (e^y), of photons from the front vicinity. This en- 
ergy is of the order of characteristic energy of 7-ray burst 
emission. 
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B. The jets of active galactic nuclei 



ACKNOWLEDGMENTS 



At radio frequencies, where VLBI can resolve the emis- 
sion regions at the miUiarcsecond scale, many of radio- 
loud active galactic nuclei (AGNs) exhibit compact jets 
(e.g., Ref. These jets are remarkably well coUi- 

mated, with opening angles typically less than a few de- 
grees. The Lorentz factor of the outflowing plasma is 
3 < To < 30. 

Recently, the Energetic Gamma Ray Experiment Tele- 
scope (EGRET) on the Compton Gamma Ray Observa- 
tory (CGRO) has detected a few ten extragalactic 
sources, which are thought to be radio galaxies favor- 
ably oriented so that the axes of the radio jets are nearly 
aligned in our observing direction (e.g., Ref. |40|). The 
main part of the 7-ray emission from the EGRET sources 
is in the range from a few ten MeV to a few GeV. Be- 
sides, very hard 7-ray s at TeV energies were detected 
from Mrk 421 with the ground-based Whipple telescope 
Hence, particles are accelerated in the jets of AGNs 
at least up to the energies of ^ 10^ — 10^^ eV. 

If the jet of AGN is free (confined solely by its own 
inertia), it resembles a conic section of a spherical wind. 
In this case, the components of the magnetic field paral- 
lel and transverse to the jet velocity are cx and 
B± cx r~^, respectively, where r is the distance from the 
central engine which is responsible for the AGN activity. 
Therefore, at large distances from the engine the field 
of the jet is mainly across its velocity, B± <C i3||, and 
our simulations may be applied to decribe the interac- 
tion between the jet front and an external medium. As 
it is shown above, in the process of such an interaction 
electrons may be accelerated to high energies, Te ::g> Fq. 
Using Eq. (p^), the mean Lorentz factor of accelerated 
electrons in the AGN frame is (Eg) ~ lO^F^. For Fq = 10, 
that is the typical value of Fq for the jets of AGNs, we 
have (Fe) - 

It was suggested |^2|-Q that the 7-ray emission of jets 
is produced when high energy electrons Compton-scatter 
an external radiation which is either the thermal UV ra- 
diation of the accretion disk around a supermassive black 
hole or the optical and UV emission of gas clouds sur- 
rounding the central engine . In both cases the typical 
energy of photons which are scattered by the high energy 
electrons is Eq '--^ 1 — 10 eV. At Fq — 10, the mean energy 
of photons after scattering is ~ eo(Fe)^ ~ 0.1 — 1 GeV, 
that is consistent with the data on 7-ray emission from 
AGNs. 



This research was supported by MINERVA Founda- 
tion, Munich / Germany. 



For Fq ^ 30, from Eq. (12) the maximum energy of 
accelerated electrons is ~ mp(?T^ ~ 10^^ eV which is 
enough to explain the energies of TeV 7-rays detected 
Ell from Mrk 421. 
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FIG. 1. Density of reversed protons in units of no in run N. 



FIG. 2. Energy spectrum of reversed protons at x < —{c/ujbp) in run N. 



FIG. 3. The beam penetration depth Xpen as a function of time t in run P (thick solid line), in run N (thin solid line), in run 
L (dotted line). 



FIG. 4. Distribution of magnetic field in run N at the moments t = (dotted line), t = 7.96Tp (thin solid line) and t = 15.9Tp 
(thick solid line). 



FIG. 5. Angular distribution of reversed protons at a; < —{c/uisp) in run N. 

FIG. 6. Maximum energy of accelerated electrons (thick line) and intensity of their synchrotron radiation per unitary area 
of the barrier front (thin line) in run N. 



FIG. 7. Energy spectrum of outflowing electrons, Vx < 0, in run N at the moment t = 15.9Tp. 



FIG. 8. Schematic drawing of a plasma beam entering a magnetic barrier. 



FIG. 9. Plots of the magnetic field Bz{x) (thin solid line), the longitudinal electric field Ex{x) (dotted line), and the x 
component of the four- velocity of electrons Ue,x {x) (thick solid line) in run SRIO. 



FIG. 10. Plots of the proton energy (solid line) and the electron energy (dotted line) in run SRIO. 

FIG. 11. Numerical resolution of the main results of the stationary case, run SRIO. o - magnetic field profile, • - electric 
field profile, x - longitudal component four-velocity of electrons. Each figure stands for five integration points. 
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FIG. 12. The value of a as a function of a/ for both the boundary condition Ex\x=xp = (thin Una) and for the boundary 
condition Ex\x=o = (thick hne). 



TABLE I. Input parameters of simulations 



Run 


Bo 


To 


a 




■^max 


^max 




(G) (103) 




(c/WBp) 


(c/uJBp) 




A 


100 


0.1 


2 


5 


10 


2.8 


B 


300 


0.3 


2 


5 


10 


2.8 


C 


300 


1 


2 


5 


12 


3.3 


D 


10* 


0.5 


2 


5 


7 


3.3 


E 


10^ 


1 


2 


10 


7 


3.5 


F 


10* 


2 


2 


5 


5 


2.4 


G 


10^ 


1 


2 


10 


5 


2.3 


H 


10* 


1 


1 


10 


5 


4.4 


I 


10^ 


1 


4 


5 


20 


3.5 


J 


10* 


1 


0.2 


5 


15 


5.1 


K 


300 


0.3 


1/3 


10 


5 


10.6 


L 


300 


0.3 


2/5 


10 


5 


10.6 


M 


300 


0.3 


4/7 


5 


30 


10.6 


N 


300 


0.3 


2/3 


5 


30 


15.9 





300 


0.3 


1 


5 


30 


10.6 


P 


300 


0.3 


2 


5 


30 


10.6 



TABLE II. Derived parameters of simulations 



Run 


(rr) 

To 


r:r»ax 

To 


ro 


(rr*) 

To 


^7 

(10-2) 


(£7) 
(MeV) 


A 


132 


603 


643 


0.61 


3.2 X 10-3 


0.013 


B 


107 


785 


809 


0.41 


0.12 


0.45 


C 


106 


464 


543 


0.57 


0.96 


5.1 


D 


119 


745 


425 


0.48 


3.4 


8.9 


E 


92 


542 


337 


0.48 


11 


14 


F 


92 


785 


235 


0.40 


12 


35 


G 


56 


402 


193 


0.43 


14 


36 


H 


116 


801 


419 


0.57 


9.9 


15 


I 


105 


427 


555 


0.41 


2.9 


9.1 


J 


9.5 


39 


41 


0.99 


0.16 


0.20 


K 


115 


573 


288 


0.94 


0.017 


0.09 


L 


424 


1473 


596 


0.73 


0.052 


0.19 


M 


296 


1038 


913 


0.69 


0.18 


0.64 


N 


489 


1275 


807 


0.50 


0.36 


0.34 





213 


1284 


821 


0.62 


0.16 


0.34 


P 


94 


550 


883 


0.63 


0.42 


0.68 
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TABLE III. Penetration of the beam particles with Fq = 300 into the magnetic barrier with Bq = 300 G for different values 
of a. 



a 2 1 2/3 4/9 2/5 1/3 
Vpen/c 0.32 0.17 0.077 0.05 



TABLE IV. Parameters of simulations of stationary collision between the beam with Fo = 300 and the magnetic barrier for 
the boundary condition Ex\x=xp = 0. 



Run 


"/ 


a 


r,, 


Te max 

Tn 


Bf 
Bn 


Ef 
B„ 


Eman 

B„ 


SRI 


0.01 


0.0098 


0.99 


1.02 


0.99 


-0.015 


-0.015 


SR2 


1/10 


0.086 


0.92 


1.2 


0.92 


-0.13 


-0.13 


SR3 


1/5 


0.15 


0.86 


1.3 


0.87 


-0.22 


-0.22 


SR4 


1/2 


0.29 


0.76 


2.1 


0.76 


-0.40 


-0.41 


SR5 


1 


0.44 


0.67 


7.1 


0.67 


-0.57 


-0.58 


SR6 


2 


0.50 


0.64 


48 


0.50 


-0.50 


-0.62 


SR7 


4 


0.51 


0.64 


97 


0.36 


-0.38 


-0.61 


SR8 


10 


0.51 


0.63 


170 


0.23 


-0.27 


-0.58 


SR9 


40 


0.51 


0.63 


276 


0.11 


-0.20 


-0.53 


SRIO 


1000 


0.51 


0.62 


434 


0.023 


-0.14 


-0.45 



TABLE V. Parameters of simulations of stationary collision between the beam with Fo = 300 and the magnetic barrier for 
the boundary condition Ex\x=o = 0. 



Run 


af 


a 


r 


Bf 
Bn 


Bn 


SLl 


0.01 


9.8 X 10"^ 


1 


0.99 


0.016 


SL2 


1/10 


0.082 


1.06 


0.90 


0.15 


SL3 


1/5 


0.13 


1.1 


0.81 


0.29 


SL4 


1/3 


0.16 


1.3 


0.69 


0.45 


SL5 


1/2 


0.15 


1.5 


0.55 


0.63 


SL6 


2/3 


0.11 


1.9 


0.41 


0.78 


SL7 


3/4 


0.085 


2.2 


0.34 


0.85 


SL8 


0.8 


0.068 


2.5 


0.29 


0.88 


SL9 


0.9 


0.032 


3.7 


0.19 


0.95 


SLIO 


1 


4.6 X 10"^ 


31 


0.02 


1 
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